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Abstract 

A new method for the determination of the real part of the elastic scattering amplitude 
is examined for high energy proton-proton at small momentum transfer. This method 
allows us to decrease the number of model assumptions, to obtain the real part in a narrow 
region of momentum transfer and to test different models. The real part is computed at 
a given point t min near t — from the known Coulomb amplitude. Hence one obtains an 
important constraint on the real part of the forward scattering amplitude and therefore 
on the p-parameter (measuring the ratio of the real to imaginary part of the scattering 
amplitude at t = 0), which can be tested at LHC. 

1 Unite de Recherche des Universites Paris 6 et Paris 7, Associee ou CNRS 



1 Introduction - The standard method for extracting 
real parts 

The standard procedure to extract the magnitude of the real part of the hadron elastic 
scattering includes a fit to the experimental data in the interference region, by minimizing 
the x 2 function: 



* 2 =£ 



(da cxp /dt(t = U)- da/dt(t = U)f 

a! 



i=l exp, 1 

where the experimental differential cross section da c * p /dt(t = ti) at the point £« and the 
statistical error A exPi ; are extracted from the measured dN/dt using, for example, the 
value of the luminosity. 

The theoretical representation of the differential cross-sections is 

^ = 2tt[|$ 1 | 2 + |$ 2 | 2 + |$ 3 | 2 + |$ 4 | 2 + 4|$ 5 | 2 ] , (2) 

where $i and $3 are the spin non-flip amplitudes. The total helicity amplitudes can be 
written as a sum of nuclear $f(s,t) and electromagnetic $f(s,t) amplitudes: 

^(s,t) = ^(s,t)+e ia ^t(t), (3) 

where $f (t) are the leading terms at high energies for the one-photon amplitude as defined, 
for example, in [T| and a is the fine-structure constant. The common phase <p is 

if = T[l + log (B(s,t)\t\/2) +v x + z/ 2 ], (4) 

where the upper (lower) sign is related to the pp (pp) scattering, B(s, t) is the slope of the 
differential cross section, 7 is the Euler constant (7 = 0.577...) and v\ and z/ 2 are small 
correction terms defining the behavior of the Coulomb-hadron phase at small momentum 
transfers (see J2] and, more recently, jH|). At very small t and fixed s, the electromagnetic 
amplitudes are such that &l(t) = $^(t) ~ a/t ,$§(£) = -®t(t) ~ a ' const. , $|(t) ~ 
—a/y/\t\ . We assume, as usual, that at high energies and small angles the one-flip and 
double-flip hadron amplitudes are small with respect to the spin-nonflip ones and that 
the hadron contributions to $1 and $3 are the same, as are the electromagnetic ones. 
Therefore 

F(s,t)=F N + F c exp(zc^). (5) 

In the 0(a) approximation, one has: 

da/dt = vr|e^F c + F N \ 2 = n [{F c + ReF N f + [a^F c + ImF N f] . (6) 

In the standard fitting procedure, one neglects the a 2 term in Eq. (JHJ) and this equation 
takes the form: 

da/dt = 7r[(F c (t)) 2 +(p(s,t) 2 + l)(ImF N (s,t)) 2 ) + 2(p(s,t) + a V (t))F c (t)ImF N (s,t)}, (7) 

where Fc(t) = =F 2aG 2 (t)/\t\ is the Coulomb amplitude (the upper sign is for pp, the lower 
sign is for pp) and G 2 (t) is the proton electromagnetic form factor squared; ReF N (s,t) 
and ImF^(s,t) are the real and imaginary parts of the hadron amplitude; p(s,t) = 
ReF]y(s, t)//mF/v(s, t). The formula ((ZJ) is used for the fit of experimental data in getting 
hadron amplitudes and the Coulomb-hadron phase in order to obtain the value of p(s, t). 



2 Computation of the real part of the spin-non-flip 
amplitude of the pp scattering from the Coulomb 
amplitude at a given point t m i n near t = 

Numerous discussions of the p-parameter (the value of p(s, t) at t = 0) measured by the 
UA4 [4] and UA4/2 [S] Collaborations in pp scattering at \fs = 541 GeV have revealed 
the ambiguity in the definition of this semi-theoretical parameter. As a result, it has been 
shown that one has some trouble in extracting from experiment the total cross sections 
and the value of the forward (t = 0) real part of the scattering amplitudes [E]-[E]- In fact, 
the problem is that we have at our disposal only one observable da/dt for two unknowns, 
the real and imaginary parts of F^(s, t). So, we need either some additional experimental 
information which would allow us to determine independently the real and imaginary parts 
of the spin non-flip hadron elastic scattering amplitude or some new ways to determine 
the magnitude of the phase of the hadron scattering amplitude with a minimum number 
of theoretical assumptions. 

Let us note two points concerning the familiar exponential forms of ReF^{s,i) and 
ImF]sf(s, t) used by experimentalists. First, for simplicity reasons, one makes the assump- 
tion that the slope of imaginary part of the scattering amplitude is equal to the slope of 
its real part in the examined range of momentum transfer, and, for the best fit, one should 
take the interval of momentum transfer sufficiently large. Second, the magnitude of p(s, t) 
thus obtained corresponds to the whole interval of momentum transfer. 

In this article, we briefly describe new and more general procedures simplifying the 
determination of elastic scattering amplitude parameters. 

From equation (jSJ, one can obtain an equation for ReF]y(s,t) for every experimental 
point v. 

ReF N (s,U) = -F C (U) ± [(l/7r)da exp /dt(U) - (a<pF c (U) + ImF N (s,ti)) 2 ] 1/2 . (8) 

The experimentalists define the imaginary part of the scattering amplitude via the 
usual exponential approximation in the small t-region 

ImF N (s, t) = a toi /(0.389 ■ 4tt) exp(£i/2), (9) 

where 0.389 is the usual converting dimensional factor for expressing <J^ Q ^ in mb. 

Equation (jHJ) shows the possibility to calculate the real part at every separate point 
ti if the imaginary part of the scattering amplitude is fixed and to check the exponential 
form of the obtained real part of the scattering amplitude (see 0). 

Let us define the sum of the real parts of the hadron and Coulomb amplitudes as 
\/Ar, so we can write: 

A£(8, U) = [ReF N (s, U) + FdU)} 2 > . (10) 

Using the differential cross sections experimental data we obtain: 

A£*(s,ti) = (I/O da exp /dt(U) - (a<pFc{U) + ImF N (s,U)) 2 . (11) 



For pp scattering at high energies, eqs. (jlllj) and (jllj) induce a remarkable property. Let 
us note that the real part of the Coulomb pp scattering amplitude is negative and exceeds 
the size of F^ p (s,t) at t — > 0, but has a large slope. As the real part of the hadron 
amplitude is known as being positive at relatively high (ISR) energies, it is obvious that 
A^ must go through zero at some value t = t^ in , i. e. 

ReF p N p (tZn) = -i^(CJ (12) 

and 

A^,CJ = 0. (13) 

Therefore A e ^ p must have a minimum at the same value t = t p £ in . 

The interpretation of Eq. (|T2j) is obvious: at fixed s, the real part of the pp amplitude 
is computed from the Coulomb amplitude at t = 1%% in - 

The magnitude of A e f^ p (s,t) as compared with A^(s, £) gives as a measure of the 
accuracy of the experiment and of the theoretical model : A e ^ p (s,t) has to be very close 
to A t ^(s, t). Consequently, A e ^ p (s,t P ^ in ) should be almost zero. The value of t^ in , defined 
in Eq. (JT3J), is determined, in fact, by the minimum of A e ^ p . If the position t p £ in of 
the minimum of A e £ p is different from the position of the zero of A^, then the model is 
questionable as concerns ReF^. This gives a powerful test for any model. 

Namely, in the case of the exponential forms, we have 



ReF p N p (s,t) 
ImF p N p (s } t) 
const = p pp (s,t min ) (14) 



/^M) = T Z£p,-J, =ff»M 



However our method gives the possibility to extract pP p {s,t^ in ) without assuming the ex- 
ponential form for ReF^^s, t), from eqs. © and (|T2j). If this numerical value of ff p (s, t pp in ) 
is significantly different from the value p pp (s, 0) extracted by a given experiment, this 
means that the exponential form of ReFff(s,t) is doubtful. 

Our method gives the possibility to extract the real part ReF^^s^i) at t = t^ in 
without assuming neither an exponential form nor any other specific form for the real 
part. Moreover, we know (e.g. from the Regge model) that the forward hadron scatter- 
ing amplitude is predominantly imaginary. Therefore a model which describes well the 
experimental dN/dt data necessarily has a good 7mFiv(s,t) for high s and small t. Our 
method precisely uses a given model for ImF^J' which is supposed to describe well the 
experimental data. 

In other words, our method is quasi model-independent : different models for ImF^{s, t) 
lead to a quite restricted range of values of t m i n . 

Let us underline, in order to avoid any misunderstandings, that our method is by 
no means aimed to extract ReFj^(s,t) as a function of t for a given s. Because of the 
dynamical dominance of the imaginary part of the hadron amplitude (described mainly by 
the Pomeron) over the real part, the expression of A e ^ p (s, ti) involves delicate cancellations 
between the two terms in Eq. (fTTj) and therefore one deals effectively with small quantities 
affected by large errors. Nothing precise about i?eF/v(s, t) as a function of t could be said 
before doing detailed and constrained fits of the data. Such fits are beyond the aim of 
the present paper. Our aim is to impose as a constraint for all existing models the zero 



in Ar(s, ti) which leads to a rather precise value of ReF N (s, t) at a special point t = t min , 
value computed from the Coulomb amplitude at the same special point. Even such a 
restricted calculation requires high-precision data and a large number of experimental 
points. The problem here is that we extract a small quantity - the real part of the hadron 
elastic amplitude - affected by large errors. In order to minimize these errors we need a 
very high-precision experiment. The only pp data we did find in literature, satisfying our 
criterium, are those at y/s = 52.8 GeV [TU]. We therefore pedagogically illustrate how 
our method works by taking the case of these data. 

In Fig. 1 we plot A^f p (s,£j) as given by eqs. (|TT|) and (JHJ), with a^ = 42.38 mb and 
B pp = 12.87 (GeV) - " ^U]- The error bars of the A e ^ p points are calculated from the 
errors bars of da exp /dt points. 

We also plot on the same figure A^(s, tj) as given by Eq. (|10p. where 

ReF p N p (s, t) = (ff p ■ (7&)/(0.389 ■ 4vr) <*p(B»ty2), (15) 

with pPP = 0.077. 

We see from Fig. 1 that there is a clear disagreement between A^(s, ti) and A e ^ p (s, ti) 
in the region 

0.03 <-t< 0.06 GeV 2 . (16) 

Namely, A t ^(s,ti) goes through zero at —t — 0.024 GeV while A e ^ p (s,ti) goes through 
a minimum at a very different value of t. Moreover, the values of the two quantities are 
very different in the region ()16j) . 

In fact the entire shape of A^ in the region (fTB|) is not consistent with the shape of 
A e ^ p . As it can be seen from Fig. 1, A^ rises very slowly in the region (|16|). while A e ^ p 
shows a rapid rise in this region 2 . 

In order to see if this discrepancy is significant we define the corresponding x 2 value: 

XIAr tr p&fm) ' [ } 

where 8 denotes the statistical error of A^f p . The overall x 2 /pt value is 2.4 for a total of 
34 points. However the major contribution to \ 2 comes from the region Ijl6|) . i.e. from 
only 10 points. Namely the value of x 2 /pt f° r f ne fi rs t 24 points is 1.2 while the value of 
the x 2 /pt f° r the last 10 points is 5.2. The effect shown in Fig. 1 is clearly statistically 
significant and can not be due to a statistical fluctuation. 

We can easily retrace the origin of the effect to da/dt itself, because of the obvious 
equality 

X \a r = X \da/dt , (18) 

which signifies that da/dt is not well fitted in the region |16j . In order to illustrate the 
effect in da/dt, we plot in Fig. 2 the quantity 

da exp Idt 
= aa i at _ 

da th /dt ' [ ' 

2 The negativity of several points of A^ p (see Fig. 1) is not important for our discussion. A very small 
correction of the normalization factor, taking into account systematical errors, and/or of the model used 
for ItuFn, eliminate this negativity. 



where we take as a theoretical model "th" the exponential model defined by eqs. (J0J) and 
(II 5|) . The quantity r is clearly different from in the region (|16j). 

In order to evaluate the Cm value we performed a polynomial fit of A^f p (s, £j) with 
the form 

a x ■ |t|" 3/2 + a 2 ■ \t\~ 2 + a 3 ■ |t| _1 + a 4 + a 5 - \t\ (\t\ - o 6 ) . (20) 

We get a x 2 /pt value of 0.73, for the following set of parameters (all parameters are 
expressed in GeV -4 ; a scale factor t = 1 GeV 2 is implicitely supposed everywhere in Eq. 
CD)): 

Oi = -0.08649, a 4 = -105.8709, 

a 2 = -0.00311, a 5 = 3154.11, (21) 

a 3 = 1.68189, a 6 = 0.04508. 

The result of this polynomial fit is shown in Fig. 3. The corresponding value of Cm is 

C„ = -0.0325 ± 0.0025 GeV 2 , (22) 

significantly different from the value t = —0.024 GeV 2 where A^(s, tj) goes through zero. 
We can therefore evaluate, from Eq. (|12jl . 



ReF^i^s = 52.8 GeV, t = CJ = 0.375 ± 0.037 (GeV) -2 (23) 

and, from Eq. (PJ), 

ImF p N p {J~s = 52.8 GeV, i = Cm) = 7.027 (GeV)" 2 . (24) 

Therefore 

pW(Vs = 52.8 GeV, t = Cm) = °- 053 ± °- 005 > ( 25 ) 

a value which is somewhat different (~ 2 standard deviations) from the value given in 
Ref. HD|: 

pPP{^~s = 52.8 GeV, t = C„) = °- 077 ± °- 009 • ( 26 ) 

The difference between the p- values, expressed by (|23j) and (|2Tj). is not highly significant, 
but it shows the power of our method in the case of high-precision experimental data. 

We verified that the influence of the specific form of the phase <fi is, as expected, small. 

The calculation presented here points out toward a real new effect revealed by our 
method. This new effect might simply mean that p is not a constant but a function of t, 
as well as B might not be a constant but also a function of t. In others words one must 
make the analysis of the experimental data with more sophisticated analytic forms of the 
scattering amplitude that the exponential one. 

The restricted range (J2~2~|) of values of t min obtained from our analysis is explicitly shown 
in Fig. 4, where we plot A^ p (y / i = 52.8 GeV, tj) computed from a model dynamically 
different from the exponential form, the Gauron-Leader-Nicolescu (GLN) model JT]. This 
model builds the scattering amplitudes from the asymptotic theorems constraints as a 
combination of Bessel functions and Regge forms, embodies the Heisenberg-Froissart In 2 s 
behavior for a? and includes the maximal Odderon ^2]- in this case, p(s,t) at a given 
s is no more a constant but varies with t. This dynamical characteristics are translated 



through the fact that A^ LAr , as it can be seen from Fig. 4, has a fast increase in the 
region (|16|). in agreement with the increase shown by A e ^ p . This fast increase shows the 
importance of ReF^ in the GLN model, in contrast with the exponential model. 

The value of t^ in , extracted from A e ^ p by using ImF^ as given by the GLN model, 
is perfectly compatible with the value (J221)- A problem still persists: the value of t^ in , 
extracted from A^, is —0.016 GeV -2 , in disagreement with the value (|22j) . 

The disagreement between A^ and A e ^ p is seen also through the values of x 2 /pt. The 
overall x 2 jpt value is comparable with the one in the exponential model case: 2.3/pt for 
a total of 34 points. Again, the major contribution comes from the last 10 points. 

It has to be noted that the GLN model has a much richer dynamical content than the 
exponential model, both from theoretical and phenomenological points of view. Moreover, 
it fits a large number of data for pp and pp scattering in a huge range of s (4.5 < \/s < 541 
GeV) values, while the exponential parameters are fixed from fits performed scattering 
by scattering and energy by energy. 

We conclude that neither the exponential model nor the GLN model can reproduce 
entirely the effect discussed in the present paper : the disagreement between A^ and 
A e ^ p . However, the stability of the value t^ in extracted from A e ^ p is remarkable: in both 
models examined in the present paper this value is perfectly compatible with the value 



There are yet not pp data at LHC. However we can evaluate A^ from Eq. (jlOj) by 
assuming an exponential form (JTKj) for ReF^, e.g. with a slope B = 22 (GeV) -2 , ot = 
111.5 mb ^3] and with p = 0.15 as illustrative values (see Fig. 5). One gets a zero in A^ 
located at t min = —0.0044 (GeV) 2 . The future small-t experiments at LHC [2] may detect 
the zero in A e ^ p leading to the computation of ReF^f at t p £ in in terms of the Coulomb 
amplitude. This would provide a strong constraint on the parameter p(^/s = 14 TeV, t = 
0). This constraint is crucial in detecting new phenomena in strong interactions (e.g., the 
Odderon presence). 

3 Conclusions 

In conclusion, we did find a new method for the determination of the real part of the 
elastic proton-proton amplitude at high s and small i at a given point t^ in near t = 0. 
The real part of the hadron amplitude is computed, at t — t^ in , from the known Coulomb 
amplitude. 

There are no hidden assumptions: we use data for dN/dt and a given form of ImFff. 
The usual method obviously needs to formulate, in addition, a given model for ReFff. 

Our method provides a powerful consistency check for the existing models and data 
and has a predictive power for the future measurements of the p-parameter at LHC. It 
requires high-precision data and a large number of experimental points. We illustrated 
how our method works by using the data at y/s = 52.8 GeV (Ref. [Toj). 

As a byproduct of our method we discovered two new effects in the data at \/s = 52.8 
GeV: 1. the significant discrepancy between A^ as defined in Eq. (IK)jl and A e ^ p as 
defined in Eq. (jlljl . A^ involving ReFjy while A e ^ p involves ImF^; 2. A e ^ p goes through 
a minimum around a t-value \t\ ~ 0.030 — 0.035 GeV -2 and has a sharp increase after 



this t- value (see Figs. 3-4). 

The dynamical origin of these general effects is still obscure. Maybe they are a result 
of oscillations in the very small t region. In order to clarify their dynamical origin, high- 
precision experimental data at a high energy other than a/s = 52.8 GeV are needed. 
In principle, the experiments which will be performed at LHC jHj could explore this 
problem. 

Stimulated by our findings, Kundrat and Lokajicek ^3] tried recently (six months 
after the publication of our results in a preprint form) to generalize our method at higher 
t-values. These authors write that the existence of a rather sharp minimum in our ap- 
proach "has provoked" them "to perform a more detailed analysis in this region with the 
help of general eikonal approach". Unfortunately, they add that our results are "bur- 
denened by two decisive discrepancies: non-allowed renormalization of experimental data 
and application of internally inconsistent simplified approach of West and Yennie" . This 
assertion is unfair, because: 1. our results are independent of any renormalization of the 
data; 2. the extension of the standard Coulomb-nuclear phase for all the range of t-values 
is beyond the scope of the present paper. Moreover, as one can see from Fig. 2 of Ref. 
[To] and from the comments of the authors on this figure, the sharp minimum in their 
generalized A R is get precisely when our equation (12) is satisfied and its locations, for 
peripheral and central behaviours, exactly correspond to our numerical results for the 
two models which we studied. In fact, the supplementary term proportional with a, in- 
duced in our equation (10) by the eikonal model of Ref. jTH] (see their equation (29)), 
produces negligible changes in the region of very small t and, therefore, our results are 
not significantly affected by the generalized formalism of Ref. [15] . 

Let us note that our method can be easily extended (with minor changes) to proton- 
antiproton scattering, by observing that, in this case, it is the combination 

ReFff - F P J (27) 

which must go through zero at some value t = t^ in . The method described in the present 
paper could be therefore used to analyze the UA4 data at \/s = 541 GeV [Sj, a complex 
work which will be done and presented in a separate paper. Of course, in general, one 
expects that t™ in ^ t™ in at fixed s. 

Our method could be also extended to the case of proton-nucleus scattering at high 
energies. 
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Figure captions 



Fig. 1 A^ (the solid curve) and A e ^ p (the triangle points) for pp scattering (Eqs. (fTUJl 
and (jlljl ) at y/s = 52.8 GeV as a function of t, computed with the exponential form of 
the amplitude (Eqs. 1J3I) and lITSI)). 



Fig. 2 The ratio r (see Eq. (J19J1 ) for pp scattering at y/s = 52.8 GeV as a function of t, 
where da th /dt is computed by using the exponential form of the amplitude (Eqs. (^J) and 



Fig. 3 Aft P (the triangle points) for pp scattering (Eq. (fTTf ) at yfs = 52.8 GeV as a 
function of £, computed with the exponential form of IitlFn (Eq. (J5|)) and fitted by the 
polynomial form (|2T?|) - (|21|) (the solid curve). The arrow indicates the position of i p r p 



mm' 



Fig. 4 A^ (the solid curve) and A e ^ p (the triangle points) for pp scattering (Eqs. (fTU|) and 
(tTTf ) at \/i = 52.8 GeV as a function oft, computed within the Gauron-Leader-Nicolescu 
(GLN) model (Ref. [Q|). 



Fig. 5 A^ for pp scattering (Eq. irU|) at yfs = 14 TeV computed within the exponential 
model (Eqs. Q and (flKjl) with the illustrative values B = 22 (GeV)" 2 , cr r = 111.5 mb 
and p = 0.15. 
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pp at Vs = 52.8 GeV 



ImF pp exponential form 
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